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Abstract 

Two- and multi-dimensional spectroscopy is used in physics and chemistry to obtain struc- 
tural and dynamical information that would otherwise be invisible by the projection into a one- 
dimensional data set such as a single emission or absorption spectrum. Here, we introduce a 
qualitatively new two-dimensional spectroscopy method by employing the carrier-envelope phase 
(CEP). Instead of measuring spectral vs. spectral information, the combined application of spec- 
tral interferometry and CEP control allows the measurement of otherwise inseparable temporal 
events on an attosecond time scale. As a specific example, we apply this general method to the 
case of attosecond pulse train generation, where it allows to separate contributions of three differ- 
ent sub-cycle electron quantum paths within one and the same laser pulse, resulting in a better 
physical understanding and quantification of the transition region between cutoff and plateau har- 
monics. The CEP-dependent separation in time between two full-cycle spaced attosecond pulses 
was determined to modulate by (54 ± 16) attoseconds. 



Two-dimensional (2D) spectroscopy is a powerful tool to study the energy-level structure, 
dynamics and couplings among quantum states in complex systems [IH3] . The measurements 
typically involve the acquisition of optical spectra as a function of dynamical (time-delay) 
coordinates in a multi-pulse nonlinear optical excitation at perturbative intensities. One of 
the time-delays is then converted into the spectral domain via Fourier transform and spec- 
tral interferometry [H [5], resulting in two-dimensional spectral vs. spectral contour plots 
that, for instance, allow separation of homogeneous from inhomogeneous spectral linewidths 
and the measurement of interstate coupling of electronic excitations. The technique was re- 
cently applied to large biomolecules, where complex energy-flow pathways could be visualized 
on an ultrafast timescale, such as the energy transfer in a photosynthetic light-harvesting 
complex j6j [7j. The first crucial step in spectral interferometry (employed in heterodyne 
2D spectroscopy) is typically the Fourier transformation of a spectral interferogram into a 
'quasi-time' domain to separate and retrieve the phase of the modulating component. So 
far, however, 2D spectroscopy methods have not been conceived to separate and analyze 
electronic events on attosecond time scales in strong few-cycle laser fields. 

In such intense few-cycle laser fields, by contrast, carrier-envelope-phase (CEP) control 
has been shown to be an important parameter, as it allows to steer electron wavepackets 
on their natural attosecond timescale [8] due to the precise temporal definition of an elec- 
tromagnetic field evolution. One problem in such experiments is that typically an electron 
wavefunction follows several quantum trajectories or transition pathways (e.g. ionizing at 
several peaks of the absolute electric field) and the experimentally observed response is a 
coherent sum over all these trajectories. Typically, electronic events separated by a half- or 
full-optical cycle interfere with each other and prevent an isolated observation of distinct 
trajectories or their quantum phases. 

Here, we present a new kind of 2D spectroscopy to quantitatively separate such intra- 
pulse strong-field processes on an attosecond time scale by combining the method of CEP 
control with spectral interferometry. This combination of techniques allows direct in-situ 
and quantitative temporal access to near-periodic strong-field electron dynamics. The re- 
sulting CEP-spectral interferometry (CEPSI) data provides additional two-dimensional tem- 
poral information in a manner similar (though qualitatively different) to the way 2D spec- 
troscopy allows the extraction of otherwise inseparable spectral features or energy-transfer 
pathways [HI [7j. Here we will apply the CEPSI method to the example where the relative 
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phase contributions of three attosecond electron trajectories (leading to an attosecond pulse 
train, APT) in a few-cycle pulse can be qualitatively and quantitatively disentangled in 
their relative timing and phase. However, the method itself can be considered to be general 
for the study of CEP-dependent interference phenomena and should thus find widespread 
use in strong-field science and physical-chemistry spectroscopy applications, especially when 
combined with an additional direct dynamical parameter such as a variable time delay of 
e.g. a pump and a probe pulse. 

APTs were previously reconstructed using cross-correlation techniques between the at- 
tosecond pulses and a moderately intense near-infrared (NIR) fundamental laser field [9j [TO] 
and theoretically developed [TTHT4] , or detected via direct optical interference of two seper- 
ately generated APTs [15J. An intrinsic property of few-cycle pulse generated APTs is the 
transistion from the adiabatic to the non-adiabatic regime [16J, where no longer all pulses in 
the train can be approximated to be equal in shape. This regime has previously been studied 
by using the non-adiabatic saddle-point approximation [T7j. More recently, a qualitative in- 
terpretation of spectral signatures due to interference of only a few attosecond pulses within 
one single few-cycle laser pulses was discussed p2J[[9], however without direct access to the 
temporal phases and relative timings which define the created APTs. Here, we illustrate 
the power of the introduced general CEPSI method by separating the phase contributions of 
individual pulses in an APT, shedding light on the underlying electron trajectory evolution. 

In the following we first employ a semi-classical model to describe the CEPSI method for 
the example of few-pulse APT generation, placing emphasis on the analysis of a non-equally 
spaced triple pulse. Subsequently we will use this model to analyze and interpret CEP- 
dependent experimental spectra. Finally, we will discuss the additional physical knowledge 
gained from our CEPSI method: it allows to connect our results to previously known regimes 
of attosecond-pulse production, while at the same time provides a better understanding of the 
spectral transition region between CEP-dependent cutoff harmonics and CEP-independent 
plateau harmonics. 

We start out with a model description of an APT: in the time domain, it can be formally 
written as a coherent sum of the electric fields of single attosecond pulses 

A (t) = A n (t ~ t n ) eX P [i$ n (t ~ t n )] , (1) 
n 

where A n (t) and $ n (t) are the amplitude and phase of each pulse, and t n denotes their 
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temporal position. In the spectral domain, interference fringes will enter the intensity distri- 
bution whenever the slow varying spectral amplitude functions A n (u) of different attosecond 
pulses share a common frequency range. Their intensity modulation can be described as 

S(lj) oc 

A n (uj)A m (uj) cos 

m>n 

where $ n (cj) is the frequency-dependent spectral phase of the n-th pulse and r nm —t m — t n 
denotes their temporal difference. A Fourier transform of this resulting spectral intensity 
modulation in Eq. ^ thus gives access to the spectral phase difference $ n (cj) — $ m (cj). 
Problems with this method arise whenever multiple contributions are overlapping, as it is 
the case for an APT with multiple half-cycle spaced (r nm = X/2) pulses. In this case, the 
contributions cannot be separated. 

In the following, we focus on the specific case of a train of three attosecond pulses gen- 
erated by (HHG) with a few-cycle strong-field driver pulse to illustrate the method. In that 
scenario, two temporal contributions exist at roughly the half-cycle spacing. The key idea 
of the CEPSI method is to separate the two overlapping temporal contributions by varying 
the CEP. As CEP variation in strong-field processes typically affects the phase of different 
contributions (e.g. multiphoton transitions [20] or here, quasi-classical electron trajectories) 
in a different way, the individual components can be separated. 

In our illustration model and in analogy to the experimental situation discussed below the 
temporal spacings Ti 2 = t 2 — ti and t 23 = £3 — £ 2 between the three pulses are chosen slightly 
asymmetrically (t 23 = Ti 2 + A), where A is on the order of the pulse duration, around 10% of 
T12. The temporal phases of the pulses are chosen to be for the first two pulses and ip for the 
third pulse. The spectral interference pattern of such a triple pulse is illustrated in Fig. [jji) 
and b). The spectra consist of alternating major and minor peaks spaced by 2n/ (ti 2 + r 23 ), 
in analogy to the diffraction pattern of a spatial triple slit. Additionally, for A ^ 0, a 
beating with frequency 2n/ A occurs. At the minima of this beating, the role of the major 
and minor peaks in the interference pattern swap which is marked by the arrows in Fig. [jji). 
This swapping occurs at different spectral positions, linearly depending on the relative phase 
ip of the third pulse with respect to the other two pulses. A Fourier analysis of the complete 
(^-dependent interference pattern both in amplitude and phase is shown in Fig.flJ). Now the 
previously defined relative pulse characteristics can be separated by observing a (^-dependent 
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FIG. 1. (color online) Spectral interference generated by a triple pulse with phases <3>i 5 2 = and 
$3 = (p and at non-equidistant temporal spacing T12 ~ I.IT23 as a function of the relative phase 
a) Spectra for different values of p. The arrows mark the spectral position where the transition 
from major to minor peaks occurs, b) 2d spectroscopic data as a continuous function of energy and 
relative phase, c) Fourier analysis of b) both in amplitude and phase where the correspondence to 
the different temporal components T12, T23 and T13 = T12 + T23 is indicated. 



amplitude and phase structure, which would not have been possible for a single value of (p 
or a (^-averaged spectrum. The amplitude peak appearing at T13 exhibits the linear phase 
difference of (p between the two pulses. The amplitude peaks at Ti 2 and r 2 3 appear very 
close to each other within their temporal bandwidth. Due to their different (^-dependence, 
a minimum appears in the temporal overlap region at (p = (2n + 1) 7r, i.e. when the relative 
phases between the first two and the latter two pulses are n out of phase. 

For the experiment, few-cycle (~ 7 fs) laser pulses at 760 nm central wavelength with 
stable CEP down to ~ 200 mrad are focused, with an / = 500 mm spherical mirror, into 
a ~ 3 mm long cell filled with neon gas at ~ 100 mbar backing pressure. To preferentially 
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phase match short trajectories, the gas cell was placed just after the laser focus [21]. The 
generated high-harmonic radiation is transmitted through a pair of 200 nm thin zirconium 
metal foils to remove direct and indirect fundamental stray light. The remaining coherent 
soft-x-ray light is spectrally analyzed with a 100 nm period free-standing Si3N4 transmission 
grating and detected with a back-illuminated x-ray-CCD camera. 

In Fig. [2^i), spatially integrated experimental spectra are shown as a function of the 
CEP. The spectra exhibit interference structures spaced by ^ 3 eV, the well known 2uj 
odd-harmonic spacing. For the energies above 95 eV the interference fringes shift linearly 
from one order to the next within one CEP cycle confirming earlier findings in the harmonic 
cutoff region [T71 [22]. Interestingly, in the energy region just below the cutoff (< 95 eV), 
a switching between different harmonic peaks occurs with a change in CEP. Repeating the 
previously illustrated Fourier analysis on these experimental data (Fig. ^p) results in a close 
agreement to the previous asymmetric three-pulse toy model, where the CEP now takes the 
role of the general phase as introduced before. A CEP dependence both for the full-cycle 
Fourier peak (713) as well as for the upper temporal part (Y23) of the half-cycle Fourier 
peak appears, whereas the lower temporal part (712) of this peak shows a constant phase as 
a function of the CEP. The local minima in the half-cylce Fourier peak in Fig. ^p) occur 
periodically at CEPs of 0.5 modulo tt. An integration over the phase for the different Fourier 
components is shown in Fig. [2^) to extract the different phase evolutions. 

To better understand the quantum-dynamical mechanism behind the formation of this 
asymmetric triple slit in the time domain, we carried out a quasi-classical strong-field ap- 
proximation simulation of HHG based on the three-step model [23j [24]. Only the short 
trajectories were considered. A 7 fs FWHM Gaussian laser pulse centered around 760 nm 
with peak electric field strength 0.11 a.u. was used in the simulations. The simulated har- 
monic spectra are shown in Fig. [3^l). At lower energies around 65 eV, the characteristic 
switching of the harmonic peaks becomes apparent, in agreement both with the experiment 
as well as with the presented initial toy model of non-equidistantly spaced triple pulses. At 
this stage it should be pointed out that propagation effects are not taken into account, which 
may be the reason behind some slight quantitative disagreement between the experiment 
and the single- atom simulations. 

Knowledge and analysis of the quasi-classical electron trajectories and their recollision en- 
ergetics allows to trace the origin of the observed triple-pulse behavior, as shown in Fig. [3b): 
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FIG. 2. (color online) Experimental results: a) Spectral interference pattern recorded as a function 
of the CEP. At energies below 95 eV a splitting of the interference structures (harmonics) occurs 
as discussed in the model of Fig. [TJ b) Fourier analysis of a): splitting of the amplitude peak as 
well as a bimodal phase behaviour is observed as a function of CEP, as in the model presented 
before, c) The phase averaged over the marked temporal regions for r\2 ( r ed triangle up), T23 (blue 
triangle down) and T13 = r\2 + T23 (circle, green) is shown to clarify the qualitatively different CEP 
dependences. 
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FIG. 3. (color online) Strong-field approximation simulation results: a) Spectral interference 

pattern as a function of the CEP. Below 75 eV, a splitting of harmonic positions at certain CEPs 

occurs, b) Fourier analysis of a) both in amplitude and phase, where only the region below 75 eV 

was considered. Again, the qualitatively different phase contributions to r\2 and T23 can be seen, 

accompanied by a minimum in the amplitude peak, c) Classically calculated electron recombination 

energies for ^cep = for the short trajectory (red solid line). The absolute value of the IR electric 

field (dashed line) and the attosecond pulse envelopes (gray shaded area) are shown. The 

different temporal spacings t\2 and T23 are indicated for the three most intense attosecond pulses 
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in the energy region, where the splitting effect occurs. 



The three most intense pulses exhibit a different temporal spacing at a given photon energy, 
indicated by the horizontal arrows Ti 2 and T23. For the two most intense pulses, a photon en- 
ergy of 65 eV corresponds to their half-cycle-plateau region (strong frequency upchirp of the 
individual attosecond pulses), whereas the third pulse contributes to this energy out of its 
half-cycle cutoff region (no frequency chirp near cutoff). The asymmetric temporal spacing 
of the three pulses thus clearly originates from the nonlinear chirp of the attosecond pulses 
that arises when the short trajectory of one pulse smoothly transits from its plateau into 
its cutoff region. It is well known [T71 [25] that the cutoff-harmonic phase depends roughly 
linearly on intensity while the plateau harmonics for the short trajectory are approximately 
independent of the intensity. By varying the CEP, we move individual attosecond pulse 
contributions up and down the envelope of the driver pulse, thus changing the intensity 
with which they are produced on a sub-cycle (attosecond) time scale. CEPSI thus allows to 
separate these different phase contributions caused by individual trajectories on a sub-cycle 
basis, surpassing previous approaches that used direct intensity variation of the entire pulse 
to extract intensity-dependent phases averaged over all half-cycles of the entire driver pulse. 
Here, the sub-cycle resolution allowed us to disentangle the trajectories in the transition 
region from the cutoff to the plateau regime within the driver pulse. 

Finally, after CEPSI allowed extracting the qualitative physical mechanism, we now 
further use it to quantify the CEP-dependent relative group delay between two full-cycle 
spaced attosecond pulses in the experiment. The full-cycle Fourier peak corresponding to 
T13 = T12 + T23 in Fig. [2)3) is fitted with a Gaussian function to determine its central position 
at each value of the CEP. The result is shown in Fig. [4| where a periodic modulation in time 
delay of (54 ± 16) attoseconds can be observed. The minima appear at CEP values where 
all three pulses are approximately equal in phase, whereas the maxima correspond to CEP 
values where phase-jumps occur across the half-cycle Fourier peak (compare to Fig. [2]). An 
explanation for this modulation is as follows: For minimal group delays, the corresponding 
CEP values belong to situations where all three pulses contribute with their plateaus. In 
contrast, when the group delay is at its maximum, the corresponding attosecond pulse con- 
figurations consist of only two pulses within their plateau region, while the third pulse on 
the trailing edge contributes with its own cutoff region (as indicated in Fig. [3)3)), resulting 
in a larger spacing T13 between pulses 2 and 3. 

In conclusion we introduced the combination of CEP sweeping and spectral interferom- 
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FIG. 4. Mean time delay T13 = T12 + T23 of the two outermost pulses in the attosecond-pulse triplet 
as a function of CEP, extracted from the experimental results in Fig. [2| The error bars denote 
the standard deviation of the Gaussian fits to the 7r-periodic data sets. A fit to the systematic 
variation in time delay results in (54 ± 16) attoseconds (horizontal dashed lines) as the third pulse 
moves into and out of the cutoff region in the spectral range considered here (70-90 eV). 

etry (CEPSI) to create a novel sub-cycle multidimensional spectroscopy method. It allows 
to separate otherwise overlapping quantum paths within the same strong-field driver pulse. 
Application of this method in HHG allows access to the transition region between the cutoff 
and plateau part in the spectrum of APTs by qualitatively and quantitatively separating 
contributions of individual sub-cycle electron trajectories. An asymmetric arrangement of 
three attosecond pulses could thus be identified and analyzed with respect to its CEP de- 
pendence and the variation of the relative timing of the two outermost pulses in the pulse 
triplet could be measured with sub-20 attosecond precision. 

While applied here to the special case of HHG as an illustration example, the method 
is general and can enhance physical insight on other strong field processes that are known 
to depend on the CEP within a single pulse, such as (above-threshold) ionization in the 
multiphoton or tunneling regime, recollision imaging, and electron localization in molecules. 
In all these processes several different quantum (transition, transport) pathways exist, with 
their phase depending sensitively on the CEP. In analogy to the high-harmonic case consid- 
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ered here, CEPSI can separate the contributing quantum paths by their phase dependence 
and thus extract novel physical mechanisms as well as quantify already existing ones at high 
temporal or spectral resolution. 
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